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Resume. We study the influence of the boundary conditions at the solid liquid 
interface on diffusion in a confined fluid. Using an hydrodynamic approach, we 
compute numerical estimates for the diffusion of a particle confined between two 
planes. Partial slip is shown to significantly influence the diffusion coefficient near 
a wall. Analytical expressions are derived in the low and high confinement limits, 
and are in good agreement with numerical results. These calculations indicate that 
diffusion of tagged particles could be used as a sensitive probe of the solid-liquid 
boundary conditions. 



1. Introduction 

The no-slip boundary condition for a fluid near a solid surface is still under 
debate [H 12]- At the macroscopic scale, the no slip boundary condition is a 
consequence of the microscopic roughness [3] . On the nanometer scale however partial 
slip is possible, and has indeed been measured experimentally [I]. This issue, which 
is important both fundamentally and for the conception of microfluidic devices, 
has motivated a number of theoretical [HI E] and numerical studies |Zj . These 
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studies have highlighted the influence of the fluid-wall interaction and pressure on 
the slippage [HI EH EE]. While chemical heterogeneities and surface roughness are 
expected to decrease slippage [TT], surfaces with special geometries can exhibit a 
" super-hydrophobic" state with a strongly increased slippage at the surface [121 EH! 
that makes fluid dynamics at solids surfaces very sensitive to surface imperfections. 
Such effects have been evidenced using micro-engineered surfaces in reference |14j . 

Slippage is usually accounted for in terms of an extrapolation length, the so- 
called slip length, here denoted as 5 [S]. This is defined as the distance inside the 
solid wall where the extrapolated flow profile vanishes. More specifically this partial 
slip boundary condition is written, for the tangential component, v t , of the velocity 
as 

with n the coordinate in the direction normal to the solid surface. The precise value 
of this slip length and its dependence on the physical and chemical characteristics 
of the surface have been investigated in a number of recent experimental studies. In 
particular, very different values for slip lengths -from a few nanometers to microns 
- have been reported using different techniques (see e.g. [I] for a review). Many of 
these techniques are indirect (pressure drop measurements [THl Ej, particle image 
velocimetry [TH], fluorescence recovery [H]), or very delicate (surface force apparatus 
|2] and Atomic force microscopy JE]). Hence the development of complementary, 
robust and non-intrusive techniques to investigate the dynamical properties of the 
solid-liquid interface would provide valuable counterparts of the previous results. 

In this manuscript, we discuss how the diffusion of tagged particles between 
walls is affected by confinement, and how such measurements could be used as a 
signature of the nature of the boundary conditions [TTI ITS} IT9]. 

We develop a theoretical and numerical approach to estimate the roles of 
confinement and slip on diffusion constrained in a planar or cylindrical pore. We make 
use of a classical hydrodynamic description, which is expected to be appropriate for 
colloidal particles, and was previously shown to be also well adapted for molecular 
diffusion 0. A numerical approach is used for the general case, and analytical 
expressions are derived in the high and low confinement regimes. 
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2. Hydro dynamic estimate of the diffusion constant 

The quantity of interest is the mean diffusion coefficient of colloidal tracers 
averaged over the measurement volume. The latter is limited here by the presence 
of the confining walls. From a theoretical point of view, a moving particle V is 
subjected to a friction force proportional to its velocity. When the motion takes 
place in a confined volume, the mobility /i depends on the boundary conditions at 
the confining walls. For a velocity U_ parallel to the boundary, the diffusion coefficient 
D\\ is given by Einstein's relation [20] 

£>U = fik B T (2) 

For a particle moving between two flat walls separated by a distance H, D\\ is a 
function of the particle radius a, the height H and the position z of the particle 
respective to the walls (in the following z will be measured by taking the origin at 
the midplane). The average diffusion coefficient in the direction parallel to the walls 
is 

The next step is to use the so called Stokes Einstein approach, i.e. to estimate 
the friction force from hydrodynamics. At low Reynolds number, the flow around the 
particle is governed by the Stokes equations : 

?/A V = VP (4) 
V.V = (5) 

where V_ is the velocity field, P the pressure field and r\ the viscosity of the fluid. 
The boundary conditions are 



fluid at rest at infinity in the unconfined directions : 




= Q 




(6) 


no slip on the particle surface : 






V\v = U 




(7) 


partial slip on solid walls, expressed by parallel V_\\ 


and perpendicular 


V±\s 


velocities and slip length 5 : 







= (8) 

V±\ s =0 (9) 

This condition is written BCs (V_) = Q. 
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The friction force experienced by the particle is then : 



F = ([ a.dS 



JJdv — 



(10) 



where a = —P I + 77 ( V U +* V U_) is the stress tensor in the fluid. 

In the next sections we provide various solutions to this boundary problem : we 
first start with a numerical "exact" solution of these equations ; the latter will be 
used subsequently as a reference solution for the approximate analytical solutions 
obtained under various assumptions. 

3. Numerical Estimates 

In this section we first start with a numerical solution of the previous equations, 
Eqs. (jlj) to (JHJ). This set of equations was solved numerically with the FEMLAB© 
software. A finite domain of size 2L x 2L x H around the particle was considered. The 
size L was chosen large enough compared to H to avoid finite size problems : typical 
values are L/a = 20 for small H, L = 3H otherwise. Space and time symmetries 
were taken into account to reduce the meshed domain for faster computations. 

From a technical point of view, the FEMLAB fluid dynamics module solves the 
bulk equations as 



with the stress tensor g_ given above. Boundary conditions are imposed according to 
Eqs. © and ©-©. 



Once the flow field has been obtained in this geometry, the force is computed 
according to equation (|IDj). Note that this way of computing the force requires a 
fine mesh since a differentiation of the velocity field is performed. A better approach 
would consist in using a weak constraint formulation so that velocity and force are 
simultaneously computed on the surface. However, such an approach is time and 
memory consuming and was not used in this work to keep computational time within 
reasonable bounds. 

Typical results are shown in figure El where the profile of the local parallel 
diffusion coefficient is plotted as a function of the altitude in the confining slab. 
As a case study, we consider the situation in which one of the two walls has a non 
zero slip length, while the no-slip boundary condition is applied at the second wall. 



V.V = 



V.ct = 



(11) 
(12) 
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Near the no slip wall, diffusion decreases from its bulk value as a result of the viscous 
friction and high velocity gradient in the fluid between the particle and the wall. This 
well known phenomenon [21] is easily explained in term of an image particle (see next 
section). For a no-slip wall, the image particle moves in the opposite direction thus 
increasing the viscous force acting on the particle. Near a partially slipping wall 
diffusion increases from the no slip case and can even be higher than the bulk value. 
In the limit 5 — > oo, diffusion reaches a high value that can be estimated using the 
image particle approach, with the image moving in the same direction as the particle 

We now turn to analytical approximate solutions of the Stokes equation in the 
previous geometry. 



4. Analytical expressions in the low confinement (large gap) limit 

When the particle is small compared to confinement height, an iterative 
reflection method can be developed, leading to an analytical expression for the friction 
force. 

In the present work we use this approach in the presence of a single, slipping, 
wall. Then, summing over the forces due to each wall yields an approximate result 
for the average diffusion coefficient. A summary of the method is given here, and 
details are discussed in Appendix El 



4-1. Reflection method with a single, slipping, wall 

The reflection method is an iterative approach [21], in which the velocity field 
V_ is expanded in the form 

V = V° + V 1 + V 2 + V 3 + ... (13) 

with each V_ n field satisfying the bulk equations (HJ)-©. The zero order field, V_q, is 
chosen as the flow field around a sphere moving in the bulk : 

Y°(r, z) = \aU (-V{z) - V(-) + £ (14) 
4 \r r 6 r 6 J 

with a the radius of the sphere. This velocity field satisfies the boundary equations 
on the particle (equation (J7J)) and at infinity (equation ©)• The method consists in 
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determining V} field such that V_° + V} satisfies the boundary conditions at infinity 
and on the solid walls ©-© : 

BC s {V l ) = -BC S (V_°) 



V 1 







Now, at this level of approximation, the boundary condition on the particle V is no 
longer satisfied by + V} and the next order Vj is defined from the reflection of 
V} on the particle as : 



V' 
V 2 



v 



V 1 



v 







The higher moments of the velocity field, V™, are built by applying iteratively the 
boundary condition on the particle and on the flat walls. 



4-2. Viscous force acting on the particle : a single wall 

The friction force experienced by the particle is the sum of individual 
contributions F n of each reflection : 

E- n = JJ dv g n -dS = JjJ^V.g n dV (15) 

where a n is the stress tensor in the fluid. For odd reflections, the velocity is regular 
in the volume of the particle. The momentum equation gives V.a n = in the domain 
occupied by V and the integral vanishes. For even reflections, the Lorentz reciprocal 
theorem [21] gives algebraically, in the limit of small particles, F n+2 = F n , 
where Vq is defined as the value of the velocity field at the center of the particle. 
One thus obtains 

F° 



u 



1 ~ u 



(16) 



with F° = 6irr]aU. As a consequence, only the velocity of the first reflected field at 
the center of the particle Vq is needed to determine mobility and diffusion coefficient. 
The calculation of this field is described in appendix A. 

Equations (|16|) and (jA.7|) give the force acting on a particle moving along a 
single planar wall as a function of the radius of the particle a, the distance from the 
wall I and the slip length S : 

Girrja 



•all 



1 - <±c 



-u 



(17) 
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C 



9 

16 (19 > 



c 



I (20) 



where the function C is defined as 

C lv) = -y - |„ - | + (J/ + ^ + §») Bfe) + ^(2v) (18) 
with -E(y) = e y Ei(y) and Ei(y) is the exponential integral function, defined as 

When 5^0 (no slip condition), one recovers the well known value 
"Z" 

derived from the method of the image particle : as mentioned above, diffusion 
decreases near a no slip wall. In the limit 5 — > oo (full slip condition), 

T 

and the presence of the wall reduces the friction force, i.e. diffusion increases, as 
measured experimentally [T7j . 

Comparisons with numerical simulations using FEMLAB© are shown in figure 
H3 Results are in good agreement down to very small distances l/a = 1.5. At large 
distance, one recovers the bulk diffusion value as expected. 

Moreover, simple and practical approximations can be obtained for the mobility 
in the limit where the distance to the wall, Z, is large compared to the slip length 5. 
Indeed an asymptotic expansion of C [y] allows to obtain 

9 1 

C [y] = — ^ (21) 

m 16i + 1 + LU 1 ; 

y Lit. 

This gives the approximate following form for the friction coefficient 

F M ^ +^£/ (22) 

1 16 1+8 

This approximation amounts to replace the distance to the wall I by Z + 5, where the 
physical meaning of the slip length in terms of an extrapolation length appears quite 
clearly in this limit. In practice, note that the expression in Eq. (|22JI leads to values 
which are within 5% to the explicit result in Eq. (fTTjl as soon as 1/8 > 0.5 ! 

After completing this work, we became aware of a similar calculation by Lauga 
and Squires [22] who computed the viscous force on a spherical particle close to a 
wall, using the same reflection method. The use of the "small particle" approximation 
corresponds to computing the flow in response to a force applied to a point-like 
particle, and can be shown to involve errors of order (a/h) 3 |2H| . 
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4-2.1. Local diffusivity in a confined geometry In order to compute the friction 
coefficient for a particle confined between two planar walls, we make the further 
assumption that each wall contributes independently to the shift in the friction force 
from its bulk value : 

E-2walls = E-lwall(. z i fi) + E.lwall{H — Z,5) — F_ bulk (23) 

where H is the distance between the two walls and here z denotes the distance to the 
bottom wall. The Einstein equation then yields for the parallel diffusion coefficient 
at a height / : 

Du = JsL ! (24) 

This expression for the friction coefficient is checked against the "exact" 
numerical results obtained using the FEMLAB software in Figures 0] and El Over 
the various slip lengths 5 and confinement gap H, the agreement is found to be quite 
good, within 6% as long as the confinement is no too strong (h/2a > 4). It can be 
observed that Eq. ()24)1 slightly underestimates the diffusion. 

Note that a different approximation could be made for the contribution of 
the two walls, by assuming that the mobility (rather than its inverse) is affected 
independently by the two walls [21] • This approximation, however, turns out to be 
less accurate than the previous approximation, in Eq. (J24I) . 

When the particle is confined to a cylindrical pore, a similar method can be 
used and provides an estimate of the viscous force acting on the particle in the low 
confinement limit (see appendix A). However, as opposed to the planar case, only 
a numerical estimation of the reflected velocity at the center of the sphere can be 
reached. An interesting difference between the planar and the cylindrical case, is that 
in the latter case, except close to the wall where the behavior is similar to a particle 
moving near a planar wall, the force acting on a tracer particle is never smaller than 
its bulk value even in the large slip length limit (i.e. the diffusion is reduced). This 
is due to the necessary recirculation of the fluid around the particle. More precisely, 
boundary conditions at infinity (no flow) imposes the overall flow rate on a section of 
the cylinder at zero. In the section centered on the particle, a negative fluid flow rate 
has to balance the positive flow rate of the particle na 2 U. Hence the viscous force 
increases from the bulk value even when 5 — > oo. In the planar geometry, recirculation 
takes place at infinity in the unconfined directions and this phenomenon does not 
take place. 
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5. Lubrication theory in the strong confinement limit 

When the confinement approaches the particle size (H ~ 2a), the main part 
of the viscous force is expected to arise from the high velocity gradient in the thin 
fluid films between each wall and the particle. In these regions, the fluid flow is 
quasi-parallel to the wall and lubrication theory [53] is expected to provide a good 
description of the velocity field. An approximation of the force acting on the particle 
can then be derived. 

We assume here that the fluid is confined between a fixed sphere and a solid wall 
moving at velocity U_ = (U, 0,0) (see figure |2J). One approximates furthermore the 
sphere by a paraboloid h(r) ~ h + |-, with r the distance to the axis of symmetry 
of the paraboloid. 

Under the lubrication assumptions [2H], the Stokes equation (jlj) for the velocity 
field, W = (W XJ W y , W x ), reduces to : 

v — ^ — = ^ll P 0> v) ( 25 ) 

The boundary conditions are written as 
dW x W x -U 



dz 5 



(26) 



dz 5 
on the wall and 



9W y _ W y 



W = (28) 
on the particle. These equations are easily integrated and using the conservation 



equation V_\\-Q,, = for the flow rate, defined as Q.. = Jq y ^ W»dz, one gets the 



following equation for the pressure, P : 

1 Mh + 46) \ (_#_ . = 

127/- 11 ^ h + 6 - 11 J ~ \2{h + S)} 
A general solution for the pressure can be written in the form P = P^ + U(r) cos(0), 
with {r,6} the angular coordinates on the planar wall. We could however not find 
an analytical solution for the previous differential equation. However a "heuristic 
solution" could be found after some manipulation of the differential equation in the 
form Il(r) = rjU rb[h(r)], with 

6 9 \n(h) 4 Hh + 6) 1 \n(h + 46) 
m ~ 5M 10~^~ + 5 o 2 + 10 P m 
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We refer to appendix IA.2I for details of the calculations leading to this result. The 
validity of this approximate expression for the pressure was checked by computing 
numerically the solution of the full differential equation (}2T?j) using a simple ODE 
solver (Mathematica©). One finds that the previous solution for the pressure differs 
from the "exact" numerical one, from less than a few percents for 5 G [0, R], 
ho E and over the full range of distance r E [0, oo[ (see figure [TT] in the 

appendix). Moreover, in the vanishing slip length limit, 5 — > 0, the previous solution 
for b[h] in Eq. (jSOjl reduces to the corresponding exact solution of the differential 
equation, which can be easily obtained as bo(h) = —-£2- 

Using this previous heuristic solution as a good approximation for the pressure, 
one may then write the force balance along the x-direction applied on the volume of 
fluid inside the cylinder r < R c (see figure |7j) as 



{f wcl ii ■+ 


/ Pn.xds) 




Jr=R c J 



At large scales R c — > 00, one may verify that the slip effect disappears : P(S,r = 
Rc) — > P(S = 0, r = Rc) and J r=Rc PdS is independent of 5 and the dependence of 
the friction force acting on the particle P come from the Rc — > 00 limit of F w all. 

A second difficulty however arises with the lubrication calculation : whatever 
the slip length 5, the friction force on the wall, F waU = J waU r]V ±W x dS, is found 
to be logarithmically divergent when R c — > 00 j2J- This can be easily verified by 
inserting in the previous friction force expression the expression W x deduced from 
the pressure field with Eq. ()30j) (see also Eq. (jB.20|) in appendix B). On the other 
hand, the difference of friction forces, AF wa u = F wa u[S] — F wa u[5 = 0], between the 
finite slip length case and the no-slip case is found to take a finite value, given in 
eq. ()B.22jl . Note that AF-p = AF wa u = AF (since the second term in Eq. (p?T|) is 
independent of 5 in the Rc — * 00 limit). 

One may however argue that this difference is a physically relevant quantity 
since the slip effects mainly affect the flow in the region with strongest confinement. 
One may indeed verify that at the lubrication level, the flows with partial slip reduces 
to the flow with the no-slip boundary condition in the region far from closest contact. 

We have plotted in Figure |H] (left) the result for AF = F wall [5] — F waU [6 = 0] 
(normalized by the bulk value of the force on the particle F^ = QirnaU) as a function 
of the minimum gap ho between the sphere and the wall. This result is compared to 
the FEMLAB calculations in the same geometry. As expected an agreement is found 
in the small gap limit, where the lubrication approximation is expected to be valid. 
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Pursuing this calculation, a diffusion coefficient can be obtained. First the 
friction coefficient on the particle situated at a distance / from the wall (with slip 
length 5) can be estimated at this level of approximation by adding to the previous 
AF the value of the friction force F wa u[l, 5 = 0} computed in the same configuration 
for a no-slip wall. Then the friction coefficient for the particle confined between two 
partially slipping wall is estimated by adding the effects of the walls on the friction 
coefficient, according to Eq. (J23J) . The mobility is finally evaluated by the inverse 
of the friction coefficient. This procedure is applied in Figure |H1 (right), where the 
numerical (FEMLAB) result has been used for the no-slip friction force F wa n[l, 5 = 0}. 
Here the diffusion coefficient is computed only for the situation where the sphere is at 
the center of the slab (z = H/2). This result is compared to "exact" results obtained 
using a full FEMLAB calculation for the particle confined between two partially 
slipping walls. Again, the lubrication approximation only yields a correct agreement 
in the small gap region, and works better for small slip lengths. When the slip length 
increases, lubrication theory overestimates the mobility : In this case, important 
contributions to the viscous force are coming from areas far from the confined zones, 
which are not properly described within the lubrication approach. 

The lubrication approach has therefore a quite limited range of application (in 
the very confined region) but the solution obtained is complementary to the low 
confinement results which works in the large gap limits. 

6. Averaged diffusivity and conclusions 

We are now in a position to compute the averaged diffusion coefficient over 
the confined slab, (D») defined in Eq. (j3J). We consider a geometry where one wall 
is characterized by a no-slip boundary condition, while a partial slip boundary 
condition, with a slip length 5 applies on the other. This configuration is chosen 
as to mimic the experimental geometry [T§] . 

Results are shown in figures El and El for various values of the confinement H 
and of the slip length 5. 

Analytical results obtained in the low confinement approximation, combined 
with the assumption of independent wall contribution, reproduce quite well the trends 
of the numerical computations. Figure H3 shows that the analytical estimate slightly 
underestimates diffusion in the low confinement limit, and tends to overestimate it 
at strong confinements. 
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In order to observe a significant dependence of diffusion on the slip length, 
two conditions are required. The particle size should not be much larger than the 
slip length, and a sufficiently strong confinement is required. With typical values of 
H ~ 4a, variations of the average diffusion constant of typically 5% to 10% would 
be expected if the slip length is changed between 0.1a and a. 

These results therefore suggest that diffusion measurements are quite sensitive 
to boundary conditions on the solid substrate. This opens new routes to measure slip 
length on the basis of the thermal motion of colloidal tracers [HI] • 



Acknowledgments L.B. thanks Yannick Almeras, with whom this work was 
initiated. 

Appendix A : First reflected field V} . 

A.l. Planar geometry 

For a particle moving at a distance I from a planar surface, the general form of 
V} satisfying equations ©-© an d © is given by : 

Vl = -L y°°y 2 " [(2 - cos 2 (w)(^| + 1)) O 1 +ikcos(u)T 1 + zk 2 cos 2 ^ 1 

x e ik(x cos(£)+° suftu))- fe| z\ ^ 1 

poo p2n j- 

Vy = — j J isin(w) \icos(u)(k\z\ + ^O 1 + -izk 2 cos^E 1 

xe ik(xcos(l)+°y S ^u))-k\z\ dkdu 



Vl = — J J [-^A;cos(n)0 1 - kT 1 +ikcos(u)(k\z\ + 1) 

x e ik(x cos{l) + °y s4u))-k\z\ dMu ^ A 3) 

where 1 , T 1 et S 1 are functions of k and u and (x,y,z) are cartesian coordinates 
centered on the particle with x along the particle velocity and z normal to the wall. 

The initial field is written in this form with (B°(/c, u) = jaU, T°(k,u) = 
— ^a 3 Uik cos(w), S°(/c, u) = 0). In the limit of a small particle, T°(/c, u) « Q°(k, u). 

Functions (6 1 , T 1 , S 1 ) are determined as the unique solution of BCs(Y° + V}) = 

: 

4 k8 + 1 
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- 3i(2lk 2 S 2 + k 2 l 2 S + kSl-S + kl 2 -l) TT . . 2kl 
= 2 2WTTT3H ^cos(.)e- 



3 (k6l + 6 + t) 



2 2OT + 1 
Integration of V^(0, 0, 0) gives 

T 
5 



3^5 



af/e 



-2jW 



7 C 



(A.5) 
(A.6) 

(A.7) 



with C\y] = -l 2 y 2 -l 2 y-l + (|y 3 + |y 2 + \y) E(y) + §y£(2y). = e^z(l, y) 
and Ei(l,y) is the exponential integral function. 



Cylindrical geometry 

For a particle moving in a cylinder, the general solution for the reflected field 
V} in (r, 0, z) cylindrical coordinates is 

' a k (X, r) cos(fc0) sin(Az) \ 



a k (X,r) 

6*(A,r) 
Cfe(A,r) 



3a£7 r+oo 
> / aA 



&A:(A, r) sin(/c0) sin(Az) 
l Cfc(A, r) cos(A;0) cos(A2;) y 



Ar 



(A.8) 



^(A)4(Ar) + *i(A)ii(Ar) + ArIli(A)i£(Ar) 

nl(\)I' k (\r) - A^i(A)J fc (Ar) - fcn£(A)/i(Ar) + AnJ(A)J fc (Ar) 



*£(A)/ fc (Ar) + Arn^(A)4(Ar) + n£(A)J fc (Ar) 



(A.9) 



where I k , I' k et I k are the first order modified Bessel functions and their derivatives. 

Bulk field ]/° for a particle moving along the cylinder axis, at a distance b from 
it, is expressed in such a form as : 

wr +~ ,+qo / r)cos(fc0)sin(Az) \ 
J/°(r,<M = ^ E ^ <*A 



«fc(A,r) 



271 fc"oo J0 



/3fe(A, r) sin(/c0) sin(Az) 
l 7fe(A, r) cos(/c0) cos(Az) y 



(A.10) 

(e r ,e ,ej 



Ar + — K k (Xr)I k (Xb) + XbK' k (Xr)I' k (Xb) 



kiK' k (Xr)I k (Xb) + -K k (Xr)I' k (Xb) 
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7fc (A,r) = 2K k (Xr)I k (Xb) + XrK' k (Xr)I k (Xb) + XbK k (Xr)I' k (Xb) 

(All) 

where -Kfc, K' k et ^ are the second order modified Bessel functions and their 
derivatives. 

f4(A), ^fc(A) and H k (X) are the unique solution of 

BC s {V° + V 1 ) = (A. 12) 

on the cylindrical wall r = R. The analytical expression is then used to compute 
numerically Vp as : 



i(7( J J . /'+ 00 

v h = -z-Y, + n *( A )) W & ) + ^n fe (A)4(Afo)] <iA (A.13) 

27r fc=-o</ 

Appendix B : Lubrication approximation 



Using reduced variables r = y / 2h ar, h = h h, 5 = h 5 and P = rill ^ 2 Jl oa p, mass 



conservation is (for compactness w remove the ~ for the reduced variables) : 

'h 3 (h + A5)„ \ „ ( h 2 



~ V„ ( ^^Y ||P + e x ..V„ ( 7T7I ^ V ) = (B.14) 



with h{r) = 1 + r 2 . 

Assuming p(r) = + rfe(r) cos(0), equation (jB.14|) becomes 
d 



Or 



d 

ra{h(r)) — (rb(r)) 
or 



+ a(h{r))b{r) + 2r(3(h{r)) = (B.15) 



with a(h) = ^g±f and m = |«. 

We could not find an exact solution for this equation. However in order to 
proceed further, we have tried to construct in a heuristic way a good approximation 
to the solution to avoid purely numerical solutions. We have proceeded as follows. 
First b is assumed to depend functionaly on h(r), as b[h{r)]. Expressing p = Poo+xblh] 
(x = r cos 6) in Eq. (jB.14|) , this equation rewrites 

- Yj| ■ (a(h) (b(h)^ + x^V^h) J + 2f3(h)x = (B.16) 

with Yii/i = {2x,2y} in cartesian coordinates. An heuristic solution is found by 
assuming a (h), /3(h) and || as constant the previous equation, which amounts to 
replace the previous equation by 
db 

-6x— + 2Ap(h)x = (B.17 
Oh 
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The constant A is adjusted so that the exact no-slip solution of Eq. (jB.14|) . 
bo(h) = — jjfs, is recovered. The solution of the Eq. (jB.17j) with A = 6/5 is 

6 9 InO) t A\n(h + 5) , 1 ln(/i + 45) 
5' 

0. Note also 



b(h) 



(B.18) 



5h5 10 5 2 5 5 2 10 5 2 

which indeed reduces to the no-slip solution bo(h) = — -gp- when 5 
that in the limit h » 5, one also recovers b(h) — > bo(h) : the pressure is independent 
of 5 far from the particle. 

The reduced viscous force acting on the wall is 



wall 



dW x 
dz 



rdrdO 



d 2 W x 



W x is determined from Stokes equation -re- 
conditions and yields 



(B.19) 

I 2 along with the boundary 



wall 



3tt 







dx 

h 2 2 1 



i(Kr) + liri(r)))— + — 



rdr + Cste(B.20) 



Deviation of the viscous force from 5 = case is, back with dimensionalized variables : 

"/i-l/, , n ,„Mh + 26)\ 2 5 1 _ 



AF = 6nr]aU 



6 



This expression can be exactly computed for the approximated b(h) given above : 

AF = ( 365 - 125 2 + IOtt V + 545 In ( \) - 545 2 In ( ~ 

360o 2 V \o J \5 

/ 1 \ 2 / 1 \ 2 

24 In (1 + 5) + 1345 In (1 + 5) 
'1 + 5" 



+ 275 2 In ( - 



- 1905 2 In (1 + 5) - 245 2 In (1 + 5f - 545 In 



5 



545^ In 



3 In (1 + 45) - 265 In (1 + 45) - 565 2 In (1 + 45) 



65 2 In (1 + 5) In (1 + 45) + 65 2 In (1 + 5) In 



1 + 45 N 
35 . 



+545 2 Li 2 



-1 



+ 65 2 Li 2 



-1-5 
35 



with Li 2 (x) the dilogarithm function defined as Li 2 (z) = J2k=i,oo z k /k 2 |2*7j . 
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1 + 5 N 
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Figure 1 : Geometry of the present calculations. A tracer particle with radius 
a diffuses in a slab with thickness H. 

Figure 2 : Numerical estimates of the reduced diffusion coefficient of a particle 
moving between a partially slipping wall (5 = 1: full line, 5 = 100 : dotted line) at 
z = and a no-slip wall at z = H, as a function of the position of the particle. From 
left to right, H/a = 3, 5, 8, 12, 17and22. 

Figure 3 : Diffusion coefficient near a single planar wall as a function of the 
distance /, for various slip length. Numerical results (solid lines) are compared with 
the analytical solution (dashed line) in the low confinement limit I » a. The slip 
length 5 increases from bottom to top. 

Figure 4 : Local diffusivity computed using the approximate analytical results 
Eq. (j24|) for 5 /a = 10 _1 (dahsed line), compared to the numerical results (solid line). 
See figure El for details and notations. 

Figure 5 : Same as in Fig. |U but for 5/a = 10 1 . See figure 121 for details and 
notations. 

Figure 6 : Flow description in the lubrication limit in the thin confined film 

Figure 7 : Sketch of the force balance in the volume r < R c . 

Figure 8 : Numerical test of the lubrication calculations : (left) plot of the 
friction force difference AF = F wa u(S) — F wa u(5 = 0) normalized by the bulk value 
Foo = QnnaU , for a single wall, as a function of the distance I to the wall. The solid 
line is the FEMLAB calculation, while the dashed line is the lubrication estimate ; 
(right) Diffusion coefficient for a particle in the middle plane of the confined geometry 
between 2 identical partially slipping walls, in the high confinement limit H/2a ~ 1. 
A good agreement between the numerical and lubrication calculations is found when 
5 -> and H/2a -> 1. 

Figure 9 : Mean diffusion coefficient between a no-slip wall and a partially 
slipping wall (5) as a function of the gap H for various slip lengths 5. Full line : 
numerical results, dashed line : low confinement approximation (average of Eq. (|24j)). 

Figure 10 : Same as figure El but plotted as a function of slip length S for 
fixed confinement H. Full line : numerical results, dashed line : low confinement 
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approximation (average of Eq. (|24Jl ). 

Figure 11 : Negative pressure — Il(r) = —rb(r) rescaled by Pq = r)U/a, as a 
function of the radial distance r. The minimum gap ho between the sphere and the 
solid surface is ho = 0.1a and the slip length is 5 = a. The solid line is the numerical 
solution of the equation for the pressure using a ODE solver (Mathematica ©). The 
dashed line is the approximate solution, Eq. ()B.18|) . see text. 
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